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Abstract. We describe transport properties of two-dimensional arrays of low 
capacitance tunnel junctions, such as the current voltage characteristic and its 
dependence on external magnetic field and temperature. We discuss several 
experiments in which the small capacitance of the junctions plays an important 
role, and we also describe the methods for fabrication and measurements. 

In arrays where the junctions have a relatively large charging energy, (i.e. 
when they have a low capacitance) and a high normal state resistance, the low 
bias resistance increases with decreasing temperature and eventually at very low 
temperature the whole array may become insulating even though the electrodes 
in the array are superconducting. This transition to the insulating state can 
be described by thermal activation, characterized by an activation energy. We 
find that for certain junction parameters the activation energy oscillates with 
magnetic field with a period corresponding to one flux quantum per unit cell. 

In an intermediate region where the junction resistance is of the order of 
the quantum resistance and the charging energy is of the order of the Joseph- 
son coupling energy, the arrays can be tuned between a superconducting and 
an insulating state with a magnetic field. We describe measurements of this 
magnetic-field-tuned superconductor insulator transition, and we show that the 
resistance data can be scaled over several orders of magnitude. Four arrays follow 
the same universal functin provided we use a modified scaling parameter. We 
find a critical exponent close to unity, in good agreement with the theory. 

At the transition the transverse (Hall) resistance is found to be very small in 
comparison with the longitudinal resistance. However, for magnetic field values 
larger than the critical value, we observe a substantial Hall resistance. The Hall 
resistance of these arrays oscillates with the applied magnetic field. Features in 
the magnetic field dependence of the Hall resistance can qualitatively be corre- 
lated to features in the derivative of the longitudinal resistance, similar to what 
is found in the quantum Hall effect. 



I INTRODUCTION 



The electrical transport in two-dimensional (2D) arrays of small Joseplison 
tunnel junctions can be described either in terms of charges or in terms of 
vortices [1,2]. Transport of charge generates a current which acts as a driving 
force for the vortices. The charge transport may be obstructed by a large 
charging energy, Ec = e^/2C, C being the capacitance of the individual junc- 
tion. On the other hand, transport of vortices generates a voltage which acts 
as a driving force for the charges. The vortex transport may be hindered if the 
Josephson coupling energy is large. At low temperature the Josephson cou- 
pling energy is given by Ej = {Rq/ R^)/ {A/2), where i?Ar is the normal state 
resistance of the individual junctions, Rq = h/Ae^ ~ 6.45 kQ is the quantum 
resistance , and A is the superconducting energy gap of the electrodes. 

Adding a single charge to an electrode in the array gives rise to an electro- 
static potential distribution which is sometimes referred to as a charge soliton 
[3]. A missing charge gives rise to the counterpart, an anti-soliton. The charge 
solitons can move in the array by single-charge tunnel events. When the elec- 
trodes are superconducting, both Single Electron Solitons (SES) and Cooper 
Pair Solitons (CPS) can exist. An exact calculation of the potential distribu- 
tion in a 49 X 47 junction array assuming only nearest neighbor coupling, is 
shown in Fig. 1, for two different cases: A single electron in the center of the 
array, giving rise to a SES, and the fundamental excitation, a soliton/anti- 



soliton pair at adjacent electrodes. The size of the soliton is given by 



as long as C > Cq. Cq is the capacitance between an electrode in the array 
and infinity. However if C is of the order of Co the nearest neighbor approach 
is not a good approximation [4]. 

On the other hand, due to the existence of a macroscopic phase of the 
superconductor, a vortex, i.e. a phase winding of 27r, with an associated 
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FIGURE 1. The exact solution for the potential distribution Uij in a 49 x 47 junction 
array, were the ratio between the junction capacitance and the self capacitance of each 

electrode is 400. The current leads arc connected at i = 1 and i = 49 a) A single electron in 
the center of the array giving rise to a single electron soliton. b) The fundamental excitation, 
a soliton/anti-soliton pair at neighboring sites. 





magnetic flux quantum $o = h/2e, can exist in a loop of the array. In a 
two-dimensional (2D) array of small Josephson tunnel junctions the duality 

between charges and vortices is especially pronounced. Therefore it is a very 
suitable object to study the dynamics of both charges and vortices. 

The arrays which have high junction resistance, Rn ^ Rq and large charg- 
ing energy, Eq ^ Ej, become insulating at low temperature due to the 
Coulomb blockade [5,6], regardless of whether the electrodes are supercon- 
ducting or normal. 

In the opposite limit, Rn ^ Rq and Ec <^ Ej, superconductivity pre- 
vails and the resistance goes to zero at low temperature. In this limit where 
the Josephson energy dominates, vortices act as classical or quasi-classical 
particles and there is a large number of papers describing such systems theo- 
retically [7-10] as well as experimentally [11-17] With normal electrodes and 
i?Af S> Rq the array goes insulating at low temperature while in the opposite 
limit Rn <^ Rq, the array stays resistive even down to very low temperature. 

In this paper we will concentrate on the arrays where the charging energy 
is comparable to or larger than the Josephson coupling energy, and where 
the junction resistance is of the order of or larger than Rq. Throughout this 
paper, we will refer to an array with normal electrodes as being in the N-state, 
and to an array with superconducting electrodes as being in the S-state. 

If the interaction between the charges (or vortices) is logarithmic, the tran- 
sition to the insulating (superconducting) state would be of the Kosterlitz- 
Thouless-Berezinskii (KTB) type [18,19]. As the temperature is increased, an 
insulating 2D array can undergo a charge unbinding transition [20] at some 
transition temperature, which results in a conductive state. Likewise, a su- 
perconducting array can undergo a vortex unbinding transition to a resistive 
state. It has been shown that the superconducting transition can be described 
as a vortex unbinding, KTB transition [11]. In more recent years there has 
also been a lot of interest in arrays where the dynamics is best described 
in terms of charges [20-27]. In several papers there has been a discussion 
whether the observed transition to the insulating state can be described as 
a KTB transition as well [20,24,25]. In a paper by Tighe et al. [25] it was 
pointed out that the transition could be well described by thermal activation 
of charge solitons. Comparing results of three different groups [20,24,25], they 
found an activation energy Ea ~ ^Ec in the N-state, a value which can be 
theoretically justified. They also suggested that in the S-state at -B = 0, 
should be jEc + Aq, where 2Ao is the superconducting energy gap at B — 
and T = 0. An interesting new result is that of Kanda and Kobayashi [27] 
where they find a thermal activation behavior at higher temperatures but a 
stronger dependence at lower temperature in the N-state. An interesting the- 
oretical development in this field is a recent paper by Feigelman et al. [28], 
where they treat the posibility of parity effects 2D arrays. In section IV we 
will describe measurements where we have investigated the transition to the 
insulating state extensively. 



For arrays in the intermediate regime {R^ ~ Rq and Ec ~ Ej) which just 
barely go superconducting, a small magnetic field can drastically change the 
low bias resistance and in fact drive the array into the insulating state [29] . In a 
theoretical description of this effect [30], the field induced excess vortices drive 
the system from a vortex glass superconducting phase into a Bose-condensed 
insulating phase. The zero-magnetic-field KTB vortex-unbinding transition is 
replaced by a field-tuned vortex-dclocalization transition. A superconductor- 
insulator(Sl) transition can also be driven by other external variables such as 
electric field [31,32] or dissipation [33]. In Section V we show experiments on 
the magnetic-field-tuned superconductor-insulator transition. The zero bias 
resistance, /?„ was measured as a function of temperature and frustration. The 
frustration, / is defined as the magnetic field normalized to Bo, where Bo is 
the field corresponding to one fiux quantum per unit cell in the array. The 
scaling curves demonstrate how Rg, plotted as a particular function of T and 
/, display a transition from insulator to superconductor. According to theory 
the resistance at the critical frustration should be universal and equal to 
-Rq. From the data of four different arrays we find a value which is of the 
order of Rq but sample dependent. We can also deduce a dynamic exponent 
close to unity, which is in agreement with the theory [30]. 

Right at the critical frustration, the Hall resistance is very small compared 
with the longitudinal resistance, indicating a small Hall effect at the SI transi- 
tion. However for frustration values larger than the critical frustration the Hall 
resistance can be relatively large. Hall measurements in both conventional su- 
perconductors [34,35] and high-Tc superconductors [36,37] have shown a sign 
reversal of the Hall resistivity in the vicinity of the superconducting tran- 
sition temperature T^, where the samples are in a mixed state. Hall mea- 
surements have also been performed on disordered superconductors near the 
superconductor-to-insulator transition [38] . In Section VI we present the frus- 
tration dependence of the longitudinal resistance Roxx and the Hall resistance 
Roxy Both the longitudinal and Hall resistance are periodic functions of the 
magnetic filed, and the Hall resistance changes sign at several magnetic fields 
within one period. We also describe the dc measurements of longitudinal 
voltage Vx and the Hall voltage as a function of bias current I^. 

II SAMPLE FABRICATION AND 
MEASUREMENTS 

The arrays are fabricated on unoxidized silicon substrates, using a com- 
bination of photo- and electron-beam-lithography and an angle evaporation 
technique. Aluminum is used for both top and bottom electrodes. The num- 
ber of junctions in each row, N, is the same as the number of junctions in each 
column, with N ranging between 10 and 168. Therefore, the array resistance 
equals the individual junction resistance, assuming a homogeneous array. 



The samples are made in two steps. First a gold contact pattern is made 
with conventional photo lithography, and then the actual array is made with 

clcctron-bcam lithography The contact pattern contains a large number of 
7x7 mm^ chips distributed over a 2 inch wafer area, each with 16 contact 
pads leading to a central area of 160x160 /xm^. A double metal layer of 20 nm 
chromium-nickel and 80 nm gold is evaporated and the redundant metal is 
lifted off in acetone. The chromium-nickel film makes the gold stick better to 
the surface. 

A double layer e-beam resist consisting of a ~210nm thick bottom layer 
of P(MMA/MAA) copolymer is spun onto the wafer and a ~60nm thick top 
layer of PMMA(950k) is used. The chip is mounted in an e-beam hthography 
instrument and the central area of each chip is exposed using the array pattern. 
A current of 20-30 pA corresponding to a beam size of about 10 nm is used. 
The beam voltage is 50 kV and the area dose is 160-200 /xC/cm^. 

Each chip is then developed in two different developers: first for ~10-20s 
in the PMMA developer which consists of a 1:3 mixture of toluene and iso- 
propanol, then for 20-40 s in the copolymer developer which consists of a 1:5 
mixture of ethyl-cellosolve-acetate (EGA) and ethanol. As an alternative the 
nontoxic mixture of 5-10 % water in isopropanol can be used to develop both 
layers in one step, with a development time of 1-2 minutes. 

After development, the resist mask contains an undercut pattern with sus- 
pended bridges [39,40] which will be used to form the junctions. By depositing 
bottom and top electrodes from different angles the overlap can be controlled. 
The base and top electrodes are evaporated from tungsten boats, while the 
substrate holder is tilted at two different angles (~ ±15 "C) to give the de- 
sired overlap. Before the top electrode is deposited, a tunnel barrier is formed 
by introducing 0.01-0.1 mbar of oxygen to the chamber for 3-10 minutes, by 
adjusting the oxidation parameters we get the desired junction resistance. 

A drawback with the angle evaporation technique is that a relatively large 
junction is formed in series with the small tunnel junction, see Fig. 3. The 
effect of this larger junction can in most cases be neglected, if its area is 
much larger than the area of the smaller junction. This is easy to make, but 
this requirement limits the minimum unit cell size of the array, which in turn 
decreases the soliton size. 

The fabrication procedure results injunctions with normal state resistances 
in the range 4 to 150 kQ, capacitances of the order of 1 fF, and Aq ~ 200 juM 
per junction. The typical unit cell size is of the order of Aceii ~ 1 A*m^ These 
values are deduced from the /^^-characteristics, assuming an offset voltage 
of Voff = Ne/2C (for a discussion of the offset value see Ref. [25]). The 
Josephson coupling energies of the individual junctions are determined [41] 
from i?jv and Aq. The superconducting transition temperature T^. for the 
aluminum is in the range of 1.35 to 1.60 K. 

The self capacitance of each electrode Cq depends on the size of the electrode 
and the size of the array, and can be estimated to be of the order 10-20 aF 



for the smaller arrays (#2 and #3) and about 2 aF for the other arrays. This 
results in a soliton size in the range of 11 to 35, measured in units of the lattice 
spacing. 

The most important parameters of the 2D arrays described in this paper 
are listed in Table 1. The arrays are divided into 3 groups. Arrays #1-3 
show a decreasing resistance for decreasing temperature, and are referred to 
as the "superconducting" arrays. They have a low resistance, Rn < Rq and 
a relatively large Josephson coupling energy, Ej/Ec > 1. 

Arrays #4-8 show also go superconducting at low temperature, but they 
display the magnetic field tuned SI transition, and are referred to as the "in- 
termediate" arrays. They have Rn ^ Rq and Ej/Ec ~ 1. 

Arrays #9-15 show an increasing resistance for decreasing temperature, 
and are referred to as the "insulating" arrays. They have a high resistance, 
Rn > 15kQ and a relatively small Josephson coupling energy, Ej/Eq < 0.5. 

The arrays are mapped onto a so called quasi-Schmid diagram [42,43] in 
Fig. 2, showing the Ej/Ec and the Rq/Rn parameters of each array. The 
Ej/Ec and Rq/Rn values in the range to oo are scaled onto the horizon- 

TABLE 1. Parameters for the 15 arrays. The resistance Rn, the capacitance C, and the 
superconducting energy gap 2Ao, were deduced from the IV-cmves. The charging energy 
Ec, and the Josephson coupling energy Ej, were calculated from these values. Bq is the 
magnetic field corresponding to one flux quantum per unit cell. A is the soliton size measured 
in units of the lattice spacing. EaN is the activation energy for normal electrodes. 
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FIGURE 2. A T = phase diagram showing the Ej/Ec and the Rq/Rn parameters for 
the measured arrays. The diagonal hne represents Aq = "^Ec- The dotted hne represents the 
border to the insulating region predicted by Fazio and Schon. The dashed line corresponds 
to a Stewart MacCumber parameter of /3c = 1. 



tal and vertical axes in the range to 2, using the function f(z)=2z/(H-z). 
The diagonal line represents Aq = 2Ec, i.e. where the charging energy for a 
Cooper-pair is equal to 2Ao. The dotted line represents the border to the in- 
sulating region predicted by Fazio and Schon [2] . The dashed line corresponds 
to a Stewart MacCumber parameter of Pc = iT^"^ /'2){Ej / Ec){Rn I RqY = 1- 
Above this line a classical Josephson junction (in the upper right part of the 
diagram) shows hysteresis [44,45]. 

Our classification of the "insulating" arrays roughly agrees with the clas- 
sification suggested by Fazio and Schon [2] shown as a dashed line in Fig. 2. 
However, a few of the arrays (#9-11) which show "insulating" behavior fall 
slightly outside the insulating area predicted by Fazio and Schon. 

In Section IV, we will concentrate on the "insulating" arrays, but we will 
also present the N-state properties of the "superconducting" arrays. The N- 
state may be thought of as lying on the x-axis of the diagram in Fig. 2. Sections 
V and VI will discuss the "intermediate" arrays. 

Fig. 3 shows a typical sample, and the inset shows the probe layout for the 
Hall measurements. There are four Hall probe pairs situated at 1/6, 1/3, 
1/2 and 5/6 the distance between the ends of the array. Each Hall probe is 



actually an array of 3 x 3 junctions, with the outer side shorted by a strip from 
which the Hall voltage is taken. In this way we can reduce the influence on the 

sample behavior arising from the presence of the voltage probes, and we are 
less sensitive to local defects possibly occurring in the probe. The resistance 
of the probe arrays is presumably similar to that of the array itself. This 
resistance is much smaller than the input impedance of our voltage amplifiers. 

The measurements are performed in a dilution refrigerator which is situated 
in an electrically shielded room. A magnetic field up to about 1400 G is applied 
perpendicular to the substrate. The magnetic field needed to produce one 
magnetic flux quantum per unit cell, Bq, varied between 10.4 G and 43.1 G 
for the different arrays. This is much less than the critical magnetic field. 
Be (~ 800 G), which was needed to bring the aluminum electrodes into the 
N-state. 

All measurements are made with the biasing and the measurement circuitry 
symmetric with respect to ground. The current- voltage(/y)-characteristics of 
the arrays are recorded at temperatures down to about 25 mK. The threshold 
voltages are deduced from the IVC measured as the voltage at which the 
current has dropped to two times the noise level (going from large bias). The 




FIGURE 3. A SEM picture of a part of an array. Each island measures 0.5 by 1 /zm, 
the overlap between neighboring islands defines the tunnel junction. The insert shows the 
probe layout. Each Hall probe consists of 3 x 3 junctions, the junctions are identical to the 
junctions in the array. 



rms current-noise determined at low voltage is typically of the order of 0.1 pA 
or less. 

For the insulating arrays we deduce the zero bias resistance Rq from the IV- 
characteristics. For the "superconducting" and "intermediate arrays" , Rq can 
also be measured with an ac-technique using lock-in amplifiers, with excitation 
currents in the range 0.01-3 nA, and frequencies of the order of 10 Hz. 

The longitudinal voltage Vx is measured at the superconducting strips at the 
ends of the array, and the Hall voltage Vy is measured on the probes located 
on opposite sites of the array. The Hall data presented here were measured 
from the probes located in the center of the array, see the inset of Fig. 3 



III CURRENT VOLTAGE CHARACTERISTICS 



The large scale /^-characteristics of the arrays are very similar for the 11 
arrays and they resemble the /F-characteristics of a single high resistance 
Josephson junction. In the S-state there is a sharp rise in the conductance at 
a voltage ~ ■ 2Ao/e. At high voltages the /^-characteristics arc linear and 
there is the usual offset voltage which is due to the Coulomb blockade [5,6]. 
In the N-state the gap feature disappears but the offset voltage remains. 

For low bias in the S-state, the /^-characteristics differ drastically for the 
different arrays. Qualitative similarities can be found between these /K-curves 
and those of single junctions biased through high impedance resistors [46]. 
Arrays #1-8 show a supercurrent-hke feature at low bias, and Rq decreases 
for decreasing temperature. Arrays #9-15 show a Coulomb blockade feature 
at low bias, and Rq increases for decreasing temperature. 



The threshold voltage Vt is the voltage at which solitons can be injected 
into the array. According to theory [3], the threshold voltage for injection of 
SESs in the N-state should be 



for a symmetrically biased array with C ^ Cq. It is important to make the 
distinction between symmetric and asymmetric (one side grounded) bias of the 
array, since the latter gives a factor of two lower threshold voltage. However, 
it has been show by Middelton and Wingreen [47] that the background charges 
modifies the picture and that the threshold voltage actually scales with the 
length of the array if the effect of random background charges is taken into 
account. 



A The threshold voltage 




(1) 



60 



J I I I I I I I I I I I I I I I I I I I I L L 




40- 



B = OG 
B= 380 G 
B = 1140 G 



20- 




VtN Vts 



-6 



-4 



' I ' ' ' I 
2 

V (mV) 



6 



FIGURE 4. The /y-characteristics of sample #15 at different values of magnetic field. 
In the S-state (B=OG), the threshold is Vts = 2.75 mV, while in the N-state {B = 1140G), 
the threshold is Vin = 0.25 mV. Note that the onset of current occurs at an even larger 
voltage for an intermediate value of B. 

In the N-state, all arrays show Coulomb blockade feature at low temperature 
and low bias, but the threshold is smeared. However there is no sharp thresh- 
old voltage in the N-state for most of the arrays. Only for array #15 , which 
has the largest Rp^, we can deduce a threshold voltage of about 0.25mV(see 
Fig. 4). Note that the onset of current occurs at a substantially higher voltage 
for intermediate magnetic fields. For those fields however the onset was more 
gradual. A similar behavior was observed in all of the "insulating" arrays. 

All of our arrays were symmetrically biased and for array #15 we get a 
theoretical value of Vt^ = 2.3 mV, according to Eq. 1. The fact that the 
measured value is substantially lower than the theoretically predicted one, is 
consistent with the picture that quantum fiuctuations effectively lower the 
energy barrier for injection of charge [26] . 

In the S-state there was a sharp threshold for all the "insulating" arrays, 
except for arrays #10 and #11. Vt is shown as a function of B for four of the 
arrays in Fig. 5. For several of the arrays Vt oscillates with B, demonstrating 
that Cooper pair solitons are injected at low magnetic fields. The period of 
oscillation corresponds to one flux quantum per unit cell and agrees well with 
the Bq values of the different arrays. The oscillations in Vt correspond to the 
oscillations in activation energy, described in Section IV, and the oscillation 
peaks in Vt and Eg, occur at the same B values (compare Figs. 8 and 5). For 




200 400 600 800 1000 200 4O0 600 800 1000 

B (G) B (G) 



FIGURE 5. The positive and negative threshold voltages as a function of B for two of 
the arrays at T » 25 mK, a) #9 At low magnetic field Vt oscillates with one flux quantum 
per unit cell, showing that Cooper pair solitons are injected, b) #15. 

increasing magnetic field the threshold voltage increases and peaks at a field 
in the range of 250 to 450 G, which is well below the critical field B^. for the 
electrodes. Vt then decreases rapidly at larger B. 

The threshold voltage in the S-state has, to our knowledge, not been de- 
scribed theoretically in the literature. However, following the arguments in 
Ref. [3], and neglecting the Josephson coupling energy, we can estimate Vt 
for two different situations. For direct injection of Cooper pair solitons, we 
get a threshold voltage which is two times higher than the value for single 
electron solitons, Vts2e = 2VtAr. The other possibility is that a Cooper-pair is 
broken up and then a single electron is injected. Then the threshold would 
be Vtse = VtN + A4Ao/e for the case of symmetric bias. For all our arrays 
Vts2e < Vtse at -B = 0, and therefore we would expect injection of Cooper 
pairs to be responsible for the threshold. 

This agrees well with our observation of the oscillating Vt at low magnetic 
field. The increase of Vt with increasing magnetic field can be understood 
in the following way. Since A, and thereby Ej, decreases with increasing 
magnetic field it becomes harder to inject Cooper-pairs, and therefore Vt in- 
creases. At some magnetic field the threshold for single electrons will become 
equal to that for Cooper-pairs {Vts2e = VtSe) and we would expect a crossover 
from Cooper-pair injection to single electron injection. This is observed in 
all the samples as a peak in Vt at magnetic fields in the range 250 to 450 G. 
Beyond this crossover VtSe decreases as a function of increasing B, due to the 
decreasing A. 

To get a more quantitative description, other effects depending on the array 
size and the background charge, as well as co-tunnehng and the Josephson 

coupling, would have to be included. The fact that the observed values at 
B = are generally lower than the "theoretical" value Vts2e, can possibly be 



explained if the Josephson coupling energy is taken into account. A step in 
this direction is a recent paper on ID-arrays where the threshold dependence 
on Ej is discussed [48]. The current above the threshold has been analyzed 
using scaling theory by Rimberg et al. [49]. 

B Hall voltages 

In an array which is strongly superconducting there is no Hall voltage since 
the Hall probes would be shorted. In an array with a strong Coulomb blockade 
the whole array is insulating and therefore the Hall probes are effectively dis- 
connected, and no Hall voltage can be measured. Therefore it is not surprising 
that it is only the "intermediate" arrays which show some Hall voltage. In 
Fig. 6a we see the IxV^ characteristics for sample #5 which shows a sharp dip 
in the current before entering the fiux-fiow regime, characterized by / being 
proportional to V. The IxVy curve, shown in Fig. 6b, displays a very similar 
feature, but the Hall voltage is about two orders of magnitude smaller. As the 
current is reversed, the Hall voltage changes sign as expected. 

In the flux-fiow regime, both and Vy increases linearly with applied cur- 
rent until the longitudinal voltage reaches the sum gap voltage of the entire 
array (A^ ■ 2Ao ~ 63 mV). At this point, the Hall voltage reaches a maximum 
of about 0.2 mV and then starts to fluctuate, gradually decreasing almost to 
zero. At still higher currents, the Hall voltage displays rich structure, which 
can also be seen in the derivative of Vx with respect to Ix- 

At very high currents, when Vx is on the normal resistance branch of the 
Ix — Vx curve, the Hall voltage Vy also increases linearly with the applied 
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FIGURE 6. The longitudinal voltage Vx (a), and the Hall voltage Vy (b) measured as a 
function of bias current for sample #5 at T « 20mK, and / « 0. Note that at low bias V^ 
and Vy are very similar 




current, with slope 2.1 Q. Because this slope is a small fraction of (0.02% 
for array #5 and 0.1% for array #7) we conclude that small non-uniformities 

exist in the array. For array ^7, the Hall voltage as a function of bias current 
shows very similar behavior to that of array #5. 

IV THE INSULATING TRANSITION 

The zero bias resistance was measured as a function of both temperature 
and magnetic field. All arrays showed an "insulating" behavior in the N-state, 
meaning that Rq increased as a function of decreasing temperature. Over a 
fairly wide range of temperature the Ro{T) dependence was exponential as can 
be seen in Fig. 7, indicating thermal activation of charge solitons. However, Rq 
saturates and is no longer temperature dependent at the lowest temperatures, 
(beyond the range displayed in Fig. 7). This saturation has been discussed in 
Ref. [26] and here we will concentrate on the thermal activation. 

In the S-state, Rq increased even more rapidly as a function of decreasing 
temperature for the "insulating" arrays. The Rq versus T curves can be di- 
vided into three temperature regions (sec Fig. 7). i) In a temperature range 
below 500 mK, Rq increases exponentially with decreasing temperature, and 
an activation energy can be defined, ii) At higher temperatures the depen- 
dence was not purely exponential due to the temperature dependence of the 
superconducting gap. iii) For the lowest temperatures Rq became larger than 
1 GQ (not shown in the figure), and could not be measured accurately with 
our experimental setup. 

In a temperature interval roughly between 200 mK and 500 mK, Rq for the 
"insulating" arrays could be fitted by a thermal activation dependence over 
the whole magnetic field range, such that 

Ro{B,T) = b-e^p^ (2) 

where Ea is the activation energy and fe is a constant. It should be noted that 
there is a region of magnetic field slightly below Be where the superconducting 

gap goes to zero in the temperature interval where the fit is made to determine 
Ea- Therefore, the Iti^Rq) vs. 1/T plots are not perfectly linear at those 
magnetic fields. However the Arrhenius law (2) is still a fair approximation. 

The fact that we do not observe a KTB charge unbinding transition in these 
arrays is not altogether surprising. It was shown by Zaikin and Panyukov [50], 
that the effect of offset charges will effectively cut of the logarithmic interaction 
between the charges. Also, the finite size of our samples is probably a limiting 
factor. 

In the N-state, Ea was found to be close to \Ec (see Table 1), and b was 

very close to i?jv, for all the "insulating" arrays. This value of Ea agrees 
well with that of Tighe et al. [25]. If the conduction is caused by thermal 



activation of single electron solitons we expect that Ea = \Ec in the N- 
state. The energy required to create a soliton anti-soliton pair starting from 
an uncharged array, i.e. the tunneling of a single charge, is the so called 



core energy Ec 



/AC = \Ec for single electrons, and four times as 



high Ecore = '^Ec for Cooper-pairs. Similarly to thermal activation in other 
systems the activation energy becomes half of the core energy. 

We next consider the S-state case, where we can imagine two alternative 
transport mechanisms. On one hand, we can calculate the energy needed to 
brake up a Cooper pair and to create a SES pair. We would expect -Ecore = 
^Ec + 2Aq and therefore, Ea = \Ec + Aq if the charge transport was entirely 
due to SESs, created from broken Cooper pairs. If on the other hand, we 
assume that only CPSs are activated, we would expect the activation energy 
to be four times higher (because of the 2e charge) than for SESs in the N-state, 
so that Ea = Ec- This picture is of course a bit naive because the Josephson 
coupling energy is not taken into account. Nevertheless, if this simple picture 
holds, we would expect that for arrays with Aq/Eq < |, thermal excitation of 
SESs would always be advantageous, and we would expect Ea = Ec + A{B). 
If on the other hand Aq/ Ec > | , we would expect a crossover from activation 
of CPSs to activation of SESs as A is suppressed by the magnetic field. All 
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FIGURE 7. The zero bias resistance Rq for array #10, vs. 1/T for B = OG, S-state 
(o), and for B = 1400 G, N-state (x). Note the very good fit to the thermal activation 
Arrhenius dependence for the N-state. The resulting slope is very close to jEc and the 
extrapolation to infinite temperature ends up right at . At low temperature the data in 
the S-state can also be fitted to an Arrhenius dependence with a slope close to jEq + Aq. 
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FIGURE 8. The activation energy Ea, as a function of magnetic field B for two of the 
arrays. At high magnetic field (in the N-state), Ea is close to jEc- The dashed lines corre- 
sponds to jEc + A(B). a) Array #13. Note the good agreement between the experimental 
data and the dashed line, b) Array #9. Note the oscillations, the period correspond to one 
flux quantum per unit cell, the amplitude is roughly equal to Ej. 



the "insulating" arrays had Aq/Eq > |, and we thus expect the crossover 
behavior with magnetic field. 

Our observations suggest a somewhat more complicated picture. As we 
enter the S-state, by lowering the magnetic field belowSc, Ea increases for all 
the "insulating" arrays. At zero magnetic field, we find that Ea is equal to 
or smaller than ^Ec + Aq but larger than Eq, for all seven arrays. In Fig. 8 
the magnetic field dependence of Ea is shown for two of the arrays, and it is 
compared to ^Ec + which is represented by a dashed line. Here, Aq 

is measured for each array and the form of the function A{B/Bc)/Ao is also 
determined from the measurements [26] . 

Prom our simple picture outlined above we would expect the ratio Aq/Eq to 
be important. We find that for arrays with a ratio Aq/ Eq less than 2, (arrays 
#12,13, and 15), Ea{B) has a behavior which is very close to lEc + A{B), 
as can be seen for array #13 in Fig. 8a. Arrays #12 and #15 show a very 
similar behavior. Tighe et al. [25] have also obtained the same result at B = 
for arrays where the ratio Aq/Ec was less than 2. However, for the other 
arrays #9-11, and #14, where Aq/ Ec was larger than 2, Ea{B) is lower than 
jEc + A{B). This can be seen in Fig. 8b, for array #9. It is obvious that the 
ratio A/Ec is important. We find that the critical value is about 2. 

In several arrays, we observe oscillations in Ea{B). The period of oscillation 
corresponds to one flux quantum per unit cell (see the inset of Fig. 8b). The 
period agrees very well with the Bq values, determined from geometry. The 
oscillations show that the Josephson couphng affects the activation energy. 
The oscillations are observed in all arrays where the measurements were taken 
with sufficiently small steps of the magnetic field. The amplitude of the oscil- 
lations is roughly Ej. Arrays with a large Ej also showed an increase of (the 



average of) Ea with increasing B resulting in a peak at 100 to 200 G (Fig. 8b). 
These effects can be understood since the creation of CPS/antiCPS pairs 

should be dependent on the Josephson coupling. For a weaker Josephson 
coupling, it should be harder to create CPS/anti-CPS pairs and Ea should 
increase. The Josephson coupling £'jcos($/$o) is affected in two ways by 
the magnetic field. At low field the cosine part is affected, resulting in an 
oscillating Ea with maxima where {B = n + ^)Bo, n being an integer. The 
oscillations of Ea demonstrate clearly that at least part of the current at low 
bias is carried by Cooper-pair solitons. At higher magnetic field the increasing 
Ea with increasing B may be explained by a decreasing A, and thereby also 
a decreasing Ej. At even higher fields A./Ec becomes smaller than 2 so that 
SES creation dominates, and thus Ea decreases with increasing field 

In summary, our results for the N-state agree well with thermal activation 
behavior, and an activation energy of ^Ec for single electron solitons can be 
extracted. For the S-state we find that as long as A/Ec < 2, pairs of SESs 
are created by breaking up Cooper pairs so that Ea = \Ec + A{B) . For larger 
values, A/Ec > 2, pairs of CPSs are responsible for a substantial part of the 
charge transport and Ea oscillates as a function of temperature. 

V THE MAGNETIC FIELD TUNED 
SUPERCONDUCTOR INSULATOR 
TRANSITION 

As mentioned previously arrays which are superconducting a low temper- 
ature can display a vortex-unbinding KTB [18] transition to a resistive state 
above a certain critical temperature T^tb- In zero magnetic field, the theory 
for KTB vortex-unbinding gives a relation between the superconducting cor- 
relation length ^ and the transition temperature Tktb- The correlation length 
is determined by a control parameter which, for example, can be the disorder 
of the system. A Josephson junction array can be described by the classi- 
cal 2D-XY model and can be associated to the KTB transition in continuous 
films [12]. In the presence of disorder, the conductivity at low temperature is 
governed by variable-range hopping(VRH) [51] and ^ should scale with Tktb- 
In a highly disordered film, the long-range vortex-pair order is destroyed, and 
Tktb is substantially suppressed compared with that of a disorder-free film. 
Provided that the transition between insulator and superconductor is contin- 
uous, right at a critical disorder, ^ should diverge and Tktb should vanish. 
If the disorder is smaller than but close to the critical disorder. Fisher [30] 
predicted, based on the analogy of VRH of vortices to VRH of electrons, that 
^ should diverge as [B — Bc)~^^ with exponent vb > 2/d = 1, where d=2 is 
the dimension of the system. The scaling theory developed for the zero field 
case implies a power law dependence of Tktb on ^-C- Tktb ~ t with 
an exponent zb of unity. Furthermore, the dual transformation suggests a 
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FIGURE 9. Linear resistance Ro at zero bias as a function of frustration for a) array 
#4 measured at several temperatures (from the top, 20, 200, 400, 600, 700 mK). At 
f = fc ~ 0.12, dRo/dT changes sign and, R* « 2.4kSl. b) array #6 at low temperatures 
(from the top 50, 70, 90, 125 mK). 

T = 0, B = Be fixed point where the magnitude of the resistivity (vector sum 
of the longitudinal and transverse components of the resistivity) should also 
be universal and equal to Rq. 

Experiments on homogeneous thin films in zero magnetic field demonstrated 
a SI transition by changing the film thickness. The transition was found to 
occur when the film sheet resistance was close to a critical value of -Rq [52]. 
A magnetic field tuned transition was found for films close to the critical 
resistance [38,53,55-57] and agreement with the scaling theory [30] was ob- 
tained. In Josephson junction arrays, a SI transition can be achieved by chang- 
ing the junction normal state resistance and the junction capacitance C 
[23,30]. The transition is found to occur near a critical point at Rj^i ~Rq 
and Ej/Ec ~ 1. We therefore expect a magnetic-field-tuned SI transition for 
arrays near the critical point. However, there is a major difference between 



TABLE 2. Parameters and scaling exponents for four of the arrays which 
showed the superconductor insulator transition. The variables are defined 
in the text. 
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Ej/Ec 


Tktb 


fc 


ZB^B 
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^oxy 






(K) 










4 


1.57 


0.70±0.04 


0.122±0.012 


4.75±0.5 


2.45±0.15 




5 


1.24 


0.46±0.04 


0.039±0.007 


8.20±1.0 


1.23±0.10 


28±5 


6 


0.80 


0.44±0.05 


0.047±0.005 


1.47±0.2 


2.20±0.15 




7 


0.90 


0.35±0.04 


0.034±0.010 


4.45±0.5 


1.61±0.15 


34±5 



a uniform film and a regular array in the presence of an external magnetic 
field. In the former case the vortices form an Abrikosov (triangular) lattice 
in the ground state [58] whereas in the latter case the vortices are pinned 
in a periodic potential, imposed by the array lattice [59]. The array lattice 
results in a ground state energy which is an oscillatory function of the applied 
magnetic field with period A/=l [59]. The ground state energy is not only 
periodic, but has minima at rational frustrations, i.e. /=l/2, 1/3, 2/3, etc. 
From one point of view, frustration is simply proportional to magnetic field 
and should be related to B in the scaling theory [30]. From another point 
of view, frustration can be consider as introducing defects from the ordered 
lattice at rational / values, in which case the KTB transition under consider- 
ation is the melting of the vortex lattice. We thus expect that the / value can 
act as a control parameter for the SI transition, and the scaling theory can 
be applied, provided one accepts that the correlation length diverges at some 
critical value of / as (/ — /c)"*^^- Experimentally, we find that the SI tran- 
sition occurs at /=n±5 and n+l/2±e (where n is an integer and 5,e « 1), 
and that the scaling analysis is applicable to both cases. 

The important parameters are listed in Table 2. The charging energies Ec 
were judged from the offset voltage of the /^-characteristics at large bias. Bg 
was about 16 G for all the intermediate arrays. For sample #6 the supercon- 
ducting mean-field-transition temperature, Tc=1.51K, was measured on an 
aluminum wire which was fabricated on the same chip as the array. The KTB 
transition temperature Tktb could be deduced from the onset of the linear 
dependence of Ro{f) [17], the values are hsted in Table 2. This method of de- 
ducing Tktb has been confirmed both theoretically [18,60] and experimentally 
[53,61]. 

Ro{f) for sample #4 at T < Tktb is depicted in Fig. 9. Below a critical 
frustration/c ~ 0.12, the resistance is lower for lower temperature, indicating 
a superconducting transition. Above fc, the resistance is higher for lower 
temperature, implying an insulating transition. The resistance at fc, R^, can 
be identified from Fig. 9a and is about 2.4 kf^. 



Fig. 9b illustrates the Ro{f) dependence for sample #6 at low temperatures. 
While arrays #4 and #5 do not show a resistance higher than their normal 
resistances at any / in the accessible temperature range, the resistance of 
arrays #6 and #7 in the most insulating case (at / ^ 0.4) are much greater 
than their R^- Remarkably, the resistance of array 7^6 changes by more than 
five orders of magnitude going from / ^ to / ~ 0.27 at 15 mK. The IV char- 
acteristics measured in the superconducting state and in the insulating state 
for array #6 are displayed in Fig. 10. In the superconducting state, the array 
exhibits clear Josephson-like current, with Ro{f ~ 0) ~ 130 fl, whereas in 
the insulating state, the array shows an insulating behavior with Ro{f ~0.27) 
>37 MQ and a back-bending feature in the IV characteristics similar to the 
behavior of higher resistance arrays {Rn ~17 kQ) reported earlier [21,23]. 
This feature is easily smeared by increasing temperature, resulting in a dras- 
tic decrease in the measured Ro in the insulating phase as seen in Fig. 9b. 
Ro{T) curves for array 7^6 in the range < / < 0.27 are shown as an inset 
of Fig. 10. The flattening-off of the resistance at non-zero frustrations at low 
temperatures is attributed to a finite size effect, explained within the context 
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FIGURE 10. Magnetic field dependence of on the IV characteristics of array #6 taken 
at 15mK. In the superconducting state (/ « 0, i?o ~ 132r2) it shows a clear supercurrent 
whereas in the insulating state (/ w 0.27, Rq «37.2MfJ) it exhibits a back-bending feature. 
The zero bias resistances differ by more than five orders of magnitude, while the change in 
magnetic field is only 4.4 G. The inset shows Ro{T) for the same array for 0< / <0.27. 



of the (vortex) VRH picture [62] . In this picture, the hopping range increases 
at low temperature. When the hopping length becomes larger than the sam- 
ple size, a temperature-independent resistance is expected. This flattening-off 
behavior was also reported in Refs. [29,26] and can be shown to depend on 
the sample size [17]. 

To appreciate the scaling theory, all the data from both the superconducting 
and the insulating sides should collapse onto a single curve when plotting 
the resistances against the scaling variable |/ — fc\ ■ j'-i/^s^'s. This is done 
using both fc and zbi^b as free parameters. The curves were determined 
by minimization of the mean square deviation from an averaged curve on 
the insulating branch. The scaling was performed in the temperature range 
50mK< T < Tktb and in the frustration range < / < 0.2. The scaling 
parameters as well as R* are listed in Table 1. Scaling curves for arrays #4 
and #6 are shown in Fig. 11. All data point in the insulating branch collapse 
onto the same trend. For the superconductor transition the low temperature 
points deviate from the general trend due to the finite size effect. 

The form of the scaling curves are very similar in the different cases, in fact 
it is possible to make the curves from all four arrays overlap by offsetting the 




FIGURE 11. i?o as a function of the scaling parameter \f — /d • j^-i/^bi^b fQj. array 
#4 and array #6, < / < 0.2. The data collapse onto one curve: the upper part for 
the insulating transition and the lower part for the superconducting transition. The low 
temperature points on the superconducting side deviate from the general trend due to the 
finite array size, a) Sample #4. The inset shows the scaling for array #4 close to full 
frustration, 0.5 < / < 0.6. b) Sample #6. The inset shows a log-log plot of the critical 
frustrations fc as a function of the Kosterlitz-Thouless transition temperature Tc for the 
four measured samples. The critical exponent zb ~1.05 can be obtained from the relation 
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FIGURE 12. The critical x- value, Xc versus the Ej/Eq value for the four arrays #4-7. 
Our data shows a clear linear dependence which is not predicted by theory 

graphs slightly in the x and the y direction. The different offsets in the y 

direction demonstrates that the zero temperature fix point is different for the 
different arrays and that wc get sample dependent values of R^. The offset in 
the x-direction shows an interesting, and unpredicted, linear dependence on 
the Ej/Ec ratio. This is shown in Fig. 12 where Xc is ploted versus the Ej/Eq 
ratio. Xc is the x-value at which the superconducting transition extrapolates 
to zero resistance. Disregarding the low temperature points on the supercon- 
ducting side, we find that we can describe all the data of the four arrays with 
a single formula (Eq. 3), where the scaling parameter deduced from theory is 
normalized to the Ej/Ec ratio. 

Ro{LT)^R;.Fi^.\I^) (3) 

Here F is a universal function describing both the insulating transition and 
the superconducting transition, and Rq is a nonuniversal, sample dependent 

parameter. 

Wc find, similar to Hebard and Paalanen [53], that the quality of the scaling 
is more sensitive to fc than to zbI'b- Using an improper /c value may cause 
large discrepancies at the left side of the plot and shghtly shift i?*, whereas an 
improper zbI'b degrades the scahng but does not change Rg* . Note that the 
scaling analysis is limited to low / values by the commensurability between 
the applied flux and the array lattice. The theory predicts that the transition 




temperature should scale with the correlation length as Tktb ~ ■ This, 
together with the fact that ^ depends on the critical frustration as fc ^ 
leads to the relation fc ~ T'^tb- Fitting fc and Tktb for the four samples 
as shown in the inset of Fig. lib, we can deduce zb = 1.05. This can be 
compared to the theoretically predicted value of unity [30]. 

In the case of /=l/2 the region of sample parameters where the phase 
transition is observable shifts to lower i?jv and larger Ej/Ec- For array #4 
we find fc = 0.5±e, with e = 0.055±0.001 and i?* is 4.33±0.15 kQ. The scaling 
is shown as an inset in Fig. 11a, the frustration range was 0.5 < / < 0.6. For 
array B, we find e = 0.030 ± 0.003 and i?* ^ 5.40 ± 0.15 kfi. Arrays #6 and 
#7 have a larger i?jv and a smaller Ej/Eq ratio, it is thus not surprising that 
the Ro(/) curves at various temperatures do not cross near / = 1/2. 

Returning to the S-1 transition at / close to zero, the scaling curves for 
all four arrays exhibit similar bifurcation shape, but with a different i?*. As 
other measurements on 2D arrays [29,54] and on superconducting films [55,56] 
also showed a different i?*, we conclude that R* is non-universal and sample 
dependent. According to theory [30] the vector sum of the longitudinal resis- 
tance Roxx (previously Ro) and the Hall resistance Roxy should be universal 
and equal to Rq at / = fc, To check this prediction, we performed measure- 
ments of Roxy for arrays #5 and #7. Four pairs of Hall voltage probes allow 
us to check the spread of the junction parameters in an array, which is found 
to be within 3%. Roxy if) shows a rich structure, the details will be published 
elsewhere. For both arrays Roxy is of the order of 30 fl, see Table 2. The 
sum of Rll^ and R^ly is thus smaller than Rqior both arrays. It should be 
noted that the smallness of i?*^^ compared with i?*^^ agrees well with recent 
experiments on thin films [38]. The Hall angle at the transition, is about 1.2° 
for both arrays. 

In the thin film case [38] a critical field was found also for the Roxy data. 
In our case the linear region in the Ix ys. Vy characteristics is very small at 
finite /, which limits the excitation current and, consequently, the resolution 
in Roxy Therefore it is hard to determine the crossing point in the Roxyif) 
curves at T < Tc for both array #5 and #7. Nevertheless, the frustration 
above which the Roxyif) curves at various T start to deviate from each other 
seems to be very close to fc- This is in contrast to the case of disordered films 
[38] where the "critical field" at which all Roxyif) curves cross is higher than 
Be and is associated with the suppression A. This is evidently not the case for 
2D arrays, since the field needed for suppression of A of our arrays is about 
800 G [26], which is much greater than the critical field (—fcBo) of a few G. 

VI THE HALL EFFECT 

Superconducting films can, to some extent, be modeled as a 2D array of 
Josephson junctions, and understanding their transport behavior can be re- 



duced to the problem of vortex dynamics. Many interesting phenomena oc- 
curring in superconducting films can also be seen in 2D Josephson junction 

arrays. In fact, the phenomena can be more easily modeled in the latter sys- 
tem because complications due to the (often unknown) microstructures do 
not exist, and phenomenological parameters such as the junction normal state 
resistance, Rjy, Josephson coupling energy Ej, and charging energy Eq, can 
be independently determined. However, there is a major difference between a 
uniform film and a regular array in the presence of an external magnetic field. 
In the former case the vortices form an Abrikosov (triangular) lattice in the 
ground state whereas in the latter case the vortices are pinned in the periodic 
potential imposed by the array lattice [59] . The array lattice of loops of area 
A, results in a ground state energy which is an oscillatory function of the 
frustration. The ground state energy is not only periodic, but has minima at 
rational frustrations, i.e. f = 1/2, 1/3,2/3, etc. [59]. In the vicinity of these 
rational frustrations, the dynamics is dominated by the motion of "defect" 
vortices [16] and the vorticity of the majority defect is reversed upon passing 
through these rational frustrations. 

The frustration dependence of Rqxx and R^y for sample #7 at various 
temperatures is shown in Fig. 13. R^x is an oscillatory function of the ap- 
plied magnetic field with period A/ = 1 and minima at rational frustrations, 
/ = 1/2, 1/3 and 2/3 as can be seen in Fig. 13. There is a critical frustration, 
fc ~ 0.034, below which the resistance decreases as the temperature is low- 
ered, and above which the resistance increases as the temperature is lowered. 
The longitudinal resistance at /c, Rq^x — 1.61 kJl , can be identified from an 
expanded view of Fig. 13a around / = fc, as well from a scaling analysis on 
these curves which we have discussed in detail in Section V, see also Ref. [63]. 

The Hall resistance Roxy as a function of frustration is shown in Fig. 13 where 
we see that it also oscillates with the applied field having the same period 
A/ =1. At / = 0, Roxy is zero, and Roxyif) has a very small negative slope. 
As the frustration is increased, Roxy goes through a minimum, increasing to 
Roxy = 34 r2 at / = fc. Thereafter it rapidly increases, reaching a maximum 
value at / ~ 0.23. As the frustration is increased further, R^xy starts to 
decrease and at / > 2/5, it becomes negative. At f ^ 0.45, Roxy reaches a local 
minimum and starts to increase to zero at / = 1/2. Roxyif) is anti-symmetric 
about / = and locally anti-symmetric about/ = 1/2. In the raw data a small 
symmetric part arises due to sample non-uniformities. This symmetric part 
can be removed by taking (Roxyif) — Roxy{—f))/'^- The removed symmetric 
part looks identical to, but is only 3% of, Rqxx- The data shown in Fig. 13 has 
the symmetric part removed. The difference in shape at / = and / = ±1 is 
not understood although it is reproducible for —1.2 < / < -f-1.2. 

The Rqxx and Roxy data can be combined to generate a third plot of the 
Hall angle © = arctan{Roxy/ Roxx) which is shown in Fig. 13 at T = 20 mK. 
Comparing our results to those of van Wees et al. [13] we find a much larger 
Hall angle. Furthermore they did not observe an anti-symmetric Roxy versus / 
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FIGURE 13. (a) The longitudinal resistance Roxx and (b) the Hall resistance Roxy, and 
(c) the Hall angle 0, as functions of frustration. Rqxx, and Roxy are shown for various 
temperatures (from the top at / «0.25, T = 20, 75, 100, 125, 150, 175mK). Rqxx is 
symmetric at / = and / = ±1/2 whereas Roxy changes sign upon passing through these 
frustrations. The data are for sample ^7 



curve. This is probably because their array was in the classical limit Ej ^ Ec-, 
whereas ours were in the quantum limit Ej ^ Ec where a finite R^xy has been 
predicted [30,64] 

The Hall angle can be interpreted as the angle between the vortex velocity 
vector Vy. and a unit vector perpendicular to the direction of current fiow J. 
The moving vortices generates an electric field given by 

E = -g„n^0oK X z (4) 

with = ±1 describing the vorticity, n„ the area density of vortices, and z 
the unit vector perpendicular to the plane of the array. Vortex motion in the 
direction of J x z creates a field E parallel to the transport current J. The 
force on a moving vortex is the Magnus force [65] 

F = Hsecpo (Vs - K) X i (5) 

where is the superfiuid electron density and K = J/n^e is the superfiuid 
velocity. Due to this force, there is a component of vortex motion parallel to 
the applied current which produces a transverse field in the —Jxz direction. 

In two dimensional regular arrays, the field E at irrational frustrations is 
generated by the motion of "defect" or vacancy vortices. From Eq. 4, it is 
clear that the sign change in the vorticity is responsible for the sign reversal 
of Roxy in the vicinity of / = 1/2, because the defects have opposite vorticity 
to the field induced vortices. 
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FIGURE 14. A comparison of Roxy (solid curve) and dRoxx/df (dotted curve) for sample 
#7 at T = 20mK, plotted with arbitrary scale. Both curves show a similar behavior, minima 
and maxima occur at the same values of frustration 



We find that the structures in Roxyif) can be correlated to structure 
dRoxx/df as shown in Fig. 14. The correlation is most clear in the vicinity 
of / = 1/2. This is similar to what is found in the Quantum Hall effect in two 
dimensional electron gases (2DEGs), where Rqxx{B) oc dRoxy/dB. That we 
find traces of the opposite derivative law is probably related to the fact that 
our system is described in termes of vortices rather than in terms of electrons 
as the QHE in 2DEGs. 

VII CONCLUSIONS 

Based on the parameters of the tunnel junction we find that we can divide 
the properties of the 2D arrays into three different categories: i) Those which 
are dominated by the Coulomb blockade and go insulating at low temperature, 
ii) those which are dominated by the Josephson effect and go superconducting 
at low temperature, and iii) the intermediate arrays which just barely become 
superconducting at low temperature, but can be made insulating by applying 
a magnetic field. 

We find that the transition for the insulating arrays can be well described 
by thermal activation of charge solitons both when the electrodes are normal 
and when they are superconducting. When the electrodes are in the normal 
state we find thermal activation of single electron solitons, with an activation 
energy fa jEc- When the electrodes are in the superconducting state we 
find a much larger Ea- For arrays with Aq/Eq < 2, the activation energy 
is simply ^Ec + A(i?) indicating that Cooper pairs are broken up and that 
pairs of single electron solitons are created. When Aq/Ec > 2, the activation 
energy oscillates with B at low magnetic field, demonstrating that Cooper 
pair solitons are created. The amplitude of these oscillations is roughly equal 
to Ej and the period corresponds to one flux quantum (h/2e) per unit cell. 

The threshold voltage for the insulating arrays also oscillates at low mag- 
netic field demonstrating that Cooper pair solitons are injected at low field. 
For increasing magnetic field the average threshold voltage increases, due to 
the decreasing Ej. In the region 250 to 450 Gauss we observe a peak in the 
threshold voltage which is interpreted as a crossover from Cooper pair soliton 
injection to single electron soliton injection. 

We have observed a frustration-tuned superconductor-insulator phase tran- 
sition in several of the "intermediate" arrays. A small applied magnetic field 
of 4.4 G can change the zero-bias resistance of an array by more than 5 orders 
of magnitude. We show scaling curves for both f=0 and /=l/2. The results 
for our four samples show a dynamic critical exponent of 1.05, in good agree- 
ment with the theory of the field-tuned S-I transition. Our data indicate a 
sample dependent R^. Moreover, we have measured the Hall resistance at fc, 
which is much smaller than Rq. 

For frustration values larger than the critical value the Hall resistance is 



substantially larger and has a rich structure as a function of applied magnetic 
field. Reversal of the sign of the Hall resistance appears at several frustrations, 

which can be attributed to the change of sign of the "defect" vortices. We 
find that the structure in Roxyif) is similar to the derivative dRoxx{f)/df. 
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